Lateral current density fronts in asymmetric double-barrier 
resonant-tunneling structures 
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We present a theoretical analysis and numerical simulations of lateral current density fronts in 
bistable resonant-tunneling diodes with Z-shaped current-voltage characteristics. The bistability is 
due to the charge accumulation in the quantum well of the double-barrier structure. We focus on 
asymmetric structures in the regime of sequential incoherent tunneling and study the dependence of 
the bistability range, the front velocity and the front width on the structure parameters. We propose 
a sectional design of a structure that is suitable for experimental observation of front propagation and 
discuss potential problems of such measurements in view of our theoretical findings. We point out 
the possibility to use sectional resonant-tunneling structures as controllable three-terminal switches. 
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I. INTRODUCTION II. MODEL OF THE DBRT STRUCTURE 
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Charge accumulation in the well of a double-barrier 
resonant tunneling (DBRT) structure leads to electro- 
static feedback that increases the energy of the quasi- 
bound state with respect to the emitter, thus making 
resonant tunneling possible for higher applied voltages ^ . 
This may result in bistability when the high current den- 
sity and low current density states coexist in a certain 
range of the applied voltage. The current-voltage char- 
acteristic becomes Z-shaped instead of the conventional 
N-shape It has been pointed out that the bistability 
provides a basis for formation of lateral current density 
patterns like traveling current density fronts . Once 
triggered, such a front leads to switching of the structure 
between the on- and off-state, depending on the applied 
voltage Front propagation in the DBRT provides 

a unique example of nonlinear pattern formation in a 
system with quantum transport. In the context of the 
DBRT studies, it represents a novel switching process 
[trigger front) which is spatially inhomogcneous and is 
likely to determine the general properties of switching 
dynamics in large-area DBRT structures. 

Recent theoretical work has essentially focused on 
the qualitative behavior of current density fronts, whereas 
the possibility of their experimental observation crucially 
depends on quantitative parameters, first of all on the 
front propagation velocity and the front width. This is 
due to the limits of temporal resolution and finite lat- 
eral dimensions of DBRT structures. In this article we 
aim to bridge this gap and present the results of numer- 
ical simulations of front dynamics in asymmetric DBRT 
structures that provide a wide range of bistability and 
thus are favorable for experimental observation. We also 
discuss the specific design of DBRT structures that are 
suitable for such experiments. 



The bistability of a DBRT structure occurs due to the 
charge accumulation in the quantum well As a result 
of the electrostatic feedback the potential of the bottom 
of the well <& depends not only on the applied voltage m, 
but also on the built-up electron concentration n (Fig.l). 
For structures with narrow well and thick barriers the 
$(w, n) dependence can be calculated by representing 
electrons in the well as a sheet charge of negligible width 
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Here -u < is the applied voltage, e < is the electron 
charge, e and ep are the relative and absolute permittiv- 
ities of the material, respectively. The effective widths 
of the emitter-well barrier bi and the well-collector bar- 
rier 62 include the half-width of the well, and 62 also in- 
cludes the width of the depleted collector layer (Fig.l) 
Eqn.(l) is also applicable in the case when the 
built-up concentration is inhomogeneously distributed in 
the transverse {x, y)— plane, provided that the transverse 
variations of n(a;, y, t) are smooth in the sense that their 
characteristic length is much larger than the effective 
thickness of the structure d. Below we will show that 
this condition is indeed met for current density fronts. 

In the regime of sequential tunneling the DBRT struc- 
ture can be described by the continuity equation for the 
electron density per unit area in the well n(a;, y, t) 3.4,6,7,9 
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where Jew{x, y) and Jwc{x, y) are local densities (per unit 
area) of the emitter-well and the well-collector currents. 
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respectively, Jj^ is the density (per unit length) of the 
transverse current in the well. 

The transverse redistribution of charge in the well is 
treated semiclassically in the drift-diffusion approxima- 
tion according to Ref. 8. The drift component of J_l is 
proportional to the lateral electrical field = — Vx$. 
Eq.(l) yields £±_ ^ V±n, and hence the drift current, 
similar to the diffusive current, turns out to be propor- 
tional to the gradient of the built-up electron concentra- 
tion. This allow us to represent the transverse current 
density in the well in terms of an effective concentration- 
dependent diffusion coefficient D{n)^: 



J_L = -eD{n)V, 
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Here /U is electron mobility in the well, rs = 
e^m) is Bohr radius in the semiconductor ma- 
terial, po = m/{TT'h?) is the 2D density of states, m is the 
effective electron mass and k is Boltzmann's constant. 
The derivation of Eq.(3) assumes a Fermi distribution in 
the well ^ . The electron temperature Tw in the well can 
be different from the substrate temperature T . In the 
following we assume a substrate temperature T = 4 K. 

The emitter-well and well-collector current densities 
are modelled by simple approximate formulas '^■^'^o 



where 
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arctan(2A/rvy) — arctan (2^1/rvK) 
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A(n, u) = E^-Ew- e$, ^(n, u) = -Ew - e$ (6) 

denote the energy of the quasibound state with respect 
to the Fermi level and the bottom of the conduction 
band in the emitter, respectively, Ef is the Fermi en- 
ergy in the emitter, Ew > E^ is the energy of the qua- 
sibound state with respect to the bottom of the well, 
Tl and are the energy broadenings associated with 
the emitter-well and the well-collector barriers, respec- 
tively, Tw = Tl + Tfl + Tscatt is the total broadening 
of the quasi-bound state resulting from the escape via 
emitter- well barrier, well-collector barrier and scattering 
in the well, respectively. The broadening of the quasi- 
bound state is characterized by a Lorentzian spectral 
function. The assumption of sequential tunneling im- 
plies that r^jF/j < Tscatt- Similar models for vertical 
transport have also been used in Refs. 6, 11. 

Eqn.(4) has a transparent physical meaning: Tl/H 
is the tunneling rate between the emitter and the well; 
[poA — n] is the difference between the number of trans- 
verse tunneling modes available for the given position 
of the quasibound state and the number of occupied 



states in the well; the last factor accoimts for smooth 
decrease of current near the resonance breaking points 
A = and = 0. Eqn.(5) describes the one-way tun- 
neling from the well to the collector with a rate Tr/Ti. 
In the stationary state the current density J is propor- 
tional to the electron concentration n stored in the well. 
Formal derivation of this transport model ^'^ is based on 
the following essential assumptions: 

(i) barriers are high in the sense that the matrix el- 
ement of the emitter-well and well-collector transitions 
do not depend on the vertical (in ^—direction) electron 
momentum, hence the barriers are fully characterized by 
constant transparencies Tl^h] 

(ii) the temperature of the substrate T and the elec- 
tron temperature of the well Tw are small compared to 
E^ and Vw 

(iii) the Fermi level in the collector is lower than the 
bottom of the well. 

Combining (l)-(5) and taking into account only one 
transverse dimension, we obtain the following nonlinear 
parabolic equation for the built-up electron charge 



/(n, u) = - {Jew{u, n) - Jwc{n)) . 
e 
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In the stationary state J = Je„ = J„c, and the local 
current-voltage characteristic J{u) = J{u,n{u)) is de- 
termined by the stationary n(u)— dependence together 
with Eqns.(4),(5). The n(u)— dependence, in turn, fol- 
lows from the steady-state condition f{n,u) = 0. The 
J{u) characteristic is bistable for sufficiently small broad- 
ening of the quasi-bound state Tw- A typical example of 
such Z-shaped characteristics is shown in Fig. 2. 

We use the following set of material and structure 
parameters: e = 12 and m = 0.067 mo (for GaAs), 
where mo is the free electron mass; E^ = 10 meV, 
^ = 10^ cmVVs, Ew = 40 meV, T = AK. The 
effective width of the emitter-well barrier is taken as 
61 = 10 nm. The effective width of the second barrier 62 
includes the width of the depletion layer in the collector 
Id ~ {eeou) / {eNjjd) (see Fig.l), which lies in the interval 
10. ..100 nm for the typical doping level Nn = 10^^ cm~^ 
in the collector 

The dependence of the transparencies Tl and on 
the actual barrier widths can be evaluated by solving 
the one-dimensional Schrodinger equation Compact 
formulas for the transparencies of GaAs/AlGaAs bar- 
riers are given in Ref. 15. Typical values are T^ fi ^ 
0.1.. .1 meV for barrier widths of the order of 10 nm 
The corresponding tunneling time is tt ~ 0.06. ..0.6 ps. 
The parameter Tscatt can be estimated on the basis of the 
electron mobility in the quantum well. Typical mobility 
values n ~ 10''...10^ cm^/Vs yield Tgcatt ~ eU/m^ 
0.1. ..1 mcV i^'i'^. The corresponding momentum re- 
laxation time Tm is of the order of 1 ps or less. Be- 
low we assume the total broadening of the quasibound 
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state to be Tw = 1---2 mcV, and take into account that 
rLjT/i < ^w- After tunneling from the emitter to the 
quantum well the transverse kinetic energy of an electron 
is transferred to the lattice. The temperature of the elec- 
tron gas in the well depends on the relation between the 
energy relaxation time Te and the tunneling time tt 
The limit cases Te <^tt and Te » tt correspond to a cold 
{Tw ~ T) and hot {Tw > T) electron gas, respectively. 
Since for resonant tunneling conditions the transverse ki- 
netic energy is limited by , the upper bound for the 
electron temperature is given by T-^ = /k « 100 K. 
Below we consider the whole range of relevant tempera- 
tures T < Tw < I k and find that the effect of electron 
heating on the front dynamics is negligible. 



III. RANGE OF BISTABILITY 

We shall start with two simple analytical estimates 
for the range of bistability Au = \uth — Uh\, where Uh 
and uth are the left and right turning points of the 
Z-shaped current-voltage characteristic that correspond 
to the holding and threshold voltages, respectively (see 
Fig. 2). In the limit Tw <^ E^ the current-voltage char- 
acteristic becomes piecewise linear. Then the ultimate 
upper bound for the bistability range Au immediately 
follows from Eqns.(4),(5) and the condition Jgyj = J^c 
(see also (A3) in the Appendix): 
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Eq.(8) indicates that low transparency Tr and large effec- 
tive thickness &2 of the well-collector barrier are favorable 
for bistability. (Note that with increasing d the intrinsic 
capacitance of the well becomes smaller, enhancing the 
electrostatic feedback from the built-up electron charge.) 

For finite broadening of the quasibound state Tw the 
bistability range shrinks or even disappears. This effect 
can be evaluated analytically by using a square-shaped 
spectral function 

Aw {E, k) = -^{e{E-Ek-Ew- e^w + Tw/2) 
i w 

-e{E-Ek-Ew-e^w-rw/2)), (9) 

where 9 is the Heaviside function, instead of the 
Lorentzian spectral function 
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that has been used to derive Eqns.(4),(5). (Here k is the 
transverse electron momentum and E)^ is the correspon- 
ding kinetic energy.) This substitution yields 
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According to (11), the effect of geometrical asymmetry 
hi 7^ 62 is not monotonic: for fixed total width d = const 
there is an optimal ratio 
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Eqn.(ll), along with Eqn.(8), predicts that Au mono- 
tonically increases with decrease of the Tr/Tl ratio. 
In contrast to this conclusion, the ranges of bistabil- 
ity calculated numerically for the full model (4), (5) (see 
Fig. 3) show that there is an optimal ratio Tr/Ti,, roughly 
Fh/Fl « 0.2, that corresponds to the widest bistability 
range. When T r is too small, the bistability vanishes due 
to the relative increase of current density in the off-state. 
Eqn.(ll) does not catch this effect because the oversim- 
plified spectral function (9) inadequately describes the 
low-branch of the current-voltage characteristic 



IV. STATIONARY PROPAGATING FRONTS 

General properties. Stationary moving fronts corre- 
spond to self-similar solutions of Eqn.(7) n{x, t) = n{x — 
vt) satisfying the boundary conditions n(— 00) = rion, 
n(-l-oo) — rioff, where rion and rioff are the electron 
concentrations in the well in the on- and off-states, re- 
spectively. Such solutions give a good approximation for 
fronts in a finite system with lateral size L much larger 
than the front width W, regardless of the actual bound- 
ary conditions for built-up concentrations at the lateral 
edges of the DBRT. A positive velocity v > corre- 
sponds to propagation of the on-state into the off-state 
(hot front), a negative velocity corresponds to propaga- 
tion of the off-state into the on-state (cold front); this 
triggers switching between the two states. In the co- 
moving frame ^ = x — vt Eqn.(7) becomes an ordinary 
differential equation 
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Multiplying by D{n)dn/d£, and integrating over /_ °° 
yields^o 
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The direction of front propagation is determined by 
the sign of the numerator in Eq.(14), and hence 
by the applied voltage u. The equal areas rule 
■^"T/ ^ ^co)D{n)dn corresponds to v = and speci- 
fies the voltage Uco at which a stationary current density 
front exists. 

Formula for the front velocity. In the limit case of small 
broadening of the quasibound state F^^ <C E^ equation 
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(13) can be integrated analytically, if we additionally as- 
sume D{n) = Do = const (see Appendix). This leads to 
the following approximate t;(u)— dependence 



TlDo {l + V^AYrB \e\{i 
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where 



A = 
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i ^ , 

OiRL = ^ = V + <^RL + 1 + VcXrl 



and Uco is determined by (A8). Generally, Eqn.(15) is 
applicable to slow fronts when |e(u — Uco)\ ^ . 

The diffusion coefficient Dq can be estimated on the 
basis of Eqn.(3), assuming drift-dominated transport in 
the wcill and choosing an average concentration (n) in the 
front wall 
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where Uoniuco) and noff(uco) are the electron concentra- 
tions in the well for on- and off-states at m = Uco, respec- 
tively. (We note that the assumption of drift-dominated 
transport is supported by the numerical results presented 
below.) Using (A4) and (A8), we substitute (16) into (15) 
and get 
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According to (8) the width of the bistability range Am 
is of order of E^ , and the second factor in (18) is of the 
order of unity for realistic structure parameters. Hence 
the order of magnitude value of the front velocity is given 

by 



respectively. The scaling rules (19) allow one to extend 
the numerical results to any value of /x and accord- 
ing to u ~ and W ~ ^J'\xjV^, provided that the 
ratio Fj^/Fi and Vw arc kept the same. 

The order of magnitude values of the front velocity 
and the front widths are 10^ cm/s (in accordance with 
(18)) and 10 /xm, respectively (Fig. 4). For lateral dimen- 
sions of the order of 100 /im, the switching times are 
thus less than 1 ns. The front width exceeds both the 
thickness of the structure d ~ 50 nm and the mean free 
path of the electron in the well which can be estimated 
as tm ~ 0.1. ..1 /im for /i = 10^ V/cms Hence the local 
approximation in (1) and the drift-diffusion approxima- 
tion (3) for the lateral transport in the well are indeed 
justified for the obtained front solutions. Cold fronts are 
faster and wider than hot fronts. This is the result of 
the higher contrast between and rio// near the right 
turning point uth of the current-voltage characteristic. 
The estimate of Eqn. (18) for the same set of param- 
eters is shown in Fig. 4(a) by the dashed line. Despite 
of the rough approximations adopted to derive (18) (see 
Appendix), it gives the correct order of the magnitude of 
the front velocity. 

In order to evaluate the significance of heating of the 
electron gas in the well we have performed simulations for 
T-w ^ Ef /k = 100 K, which corresponds to the upper 
bound of the electron temperature Tw (Fig. 4a, curve 2). 
It is found that heating has a visible but minor effect on 
front velocity. This indicates the drift-dominated lateral 
transport in the quantum well. To confirm this conclu- 
sion we have also performed simulations for = 0, 
when lateral diffusion in the well is minimized and have 
found that the front velocity is practically the same as for 
Tw = T = 4K (Fig. 4a, curve 1). Generally, according 
to (3) the drift component of the lateral current is much 
larger than the diffusion component for n » PokTw ^■ 
As it follows from the characteristic shown in Fig. 2, for 
T = Tw this condition is satisfied for almost the whole 



range of concentrations [n 



on : ^of f 



within the propagat- 



ing front. Note that since the effective diffusion coeffi- 
cient in (3) is concentration-dependent, the voltage Uco 
shifts with electron temperature Tw- The right and left 
extreme points for Uco correspond to the limits Tw — > 
(pure drift) and Tw oo (pure diffusion), respectively. 



V ~ 
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(18) 



We get V ~ 10^ cm/s for E^ = 10 meV, F^, = 1 meV and 
H = 10^ cmVVs. 

Numerical results. Numerical results for the front ve- 
locity v{u) and front width W{u) are summarized in 
Fig. 4. It follows from Eqn. (7) that these quantities scale 



as 



^/Df and 




(19) 



DISCUSSION: TOWARDS EXPERIMENTAL 
OBSERVATION OF FRONT DYNAMICS. 



Since the pioneering work ^ the experimental obser- 
vation of bistability in DBRT structures has become a 
routine. Despite of this, to the best of our knowledge 
lateral current density patterns in the DBRT, predicted 
already several years ago, have not yet been observed ex- 
perimentally. Below we discuss the potential problems of 
such measurements in view of our theoretical findings. 



4 



Sectional design of DBRT strueture. Achievements of 
modern growth technology have made it possible to man- 
ufacture resonant-tunneling structures with sophisticated 
design 21.22^ opening up the way to study transverse cur- 
rent density patterns in the DBRT experimentally. Here 
we propose the simplest design of a DBRT diode with a 
sectional emitter that allows to trigger a current density 
front and to observe its propagation by means of electri- 
cal measurements. Such a device constitutes of identical 
sections with individual emitters but common collector 
so that the bias can be applied separately to each section. 
To trigger a hot front the system should be prepared in 
the uniform off-state. This can be done by applying a 
voltage u > Uh to all emitter sections and subsequently 
decreasing this voltage to a value Uh < u < Uco- Then 
the DBRT can be switched on locally by a short negative 
voltage pulse u < Uh, applied to either lateral or central 
section of the device. The resulting front propagation 
can be monitored by measuring the electrical current in 
the parallel sections of the device. 

Barriers asymmetry. Our results suggest that to 
achieve the widest bistability range the transparency of 
the well-collector barrier should be approximately 5 times 
smaller than the transparency of the emitter- well barrier. 

Front velocity and time resolution. Our simulations 
show that the order of magnitude values of the front ve- 
locity and width are lO'^cm/s and 10 ^m, respectively. 
It implies that the lateral dimension of suitable DBRT 
structure should exceed 100 /xm. Since it takes about 1 
ns for the front to cross such a structure, reliable mea- 
surements of front propagation demand for time resolu- 
tion of 100 ps or better. The fronts become slower with 
decreasing transparency of the barriers Tl and Tr. How- 
ever, then the front width W will increase (see Eq.(19)), 
setting severe additional limitations on the lateral dimen- 
sion of the appropriate DBRT structure. For example, 
decreasing Tl^r by one order of magnitude to the values 
Tl = 0.05 meV, = 0.01 meV results in a front width of 
the order of ^ 30 /im, whereas the width of the cold front 
can exceed 100 /xm. Finally we note that in the regime of 
coherent tunneling the front velocity is expected to be of 
the order of the Fermi velocity v j = \/'2E^ /m* ^ . This 
yields ~ 10^ cm/s for = 10 meV. 

Hot or cold fronts ? The width of a hot front is con- 
siderably smaller than the width of cold fronts and its 
velocity is somewhat lower (Fig. 4). Hence a hot front is 
expected to be an easier target for experimental observa- 
tion than a cold front. 

Pinning of the front. Pinning on small imperfections 
of the reference media has been recognized as an impor- 
tant mechanism capable of preventing the front propa- 
gation in nonlinear systems In the case under con- 
sideration both technological imperfections of the DBRT 
structure and the spatial inhomogeneity due to the sec- 
tional DBRT design may play a role. Due to the exponen- 
tial dependence of the barrier transparencies F^.j? on the 
barrier width, the embedded fluctuations of the barrier 
width are the most important among technological im- 



perfections. In presence of imperfections the stationary 
current-voltage characteristic includes additional stable 
branches that correspond to nonuniform current density 
profiles associated with imperfections These steady 
non-uniform states remain observable, though the uni- 
formly propagating fronts may not be possible in low- 
quality structures. 

Effect of the external circuit. Generally, bistable 
DBRT elements can be operated via an external load. In 
contrast to the well known case of S-type bistability, in 
the system with Z-type bistability an external load pro- 
vides a positive, not negative feedback upon front dynam- 
ics. In particular, stationary current density patterns are 
not stable even under current-controlled conditions ^. An 
active external circuit which simulates a negative exter- 
nal load resistance and reverses the sign of feedback has 
been implemented experimentally in Rcf. 2. We refer to 
Refs. 7, 25 for the further discussion of the effect of exter- 
nal circuit on dynamics and stability of current density 
patterns. 

DBRT-based 3-electrode switch. The sectional DBRT 
with two small (control) sections separated by a large 
(main) section can be potentially used in applications 
as a controllable switch. By applying a control voltage 
pulse to the side sections, the rest of the structure can 
be switched to the on- or off-state via propagation of a 
lateral current density front. Note that in semiconductor 
device physics propagation of switching fronts in three- 
electrode bistable devices is a well known process which 
is used, e.g, in thyristors in certain regimes of operation 

26,27 
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APPENDIX: ANALYTICAL APPROXIMATION 
FOR THE FRONT VELOCITY. 



In the limit Tw <C Ep the local kinetic function be- 



comes piece wise linear: 
f{n, u) = —(A + aRL + 1) 



for n e 
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where the following notations have been introduced 
46162 



A = 



— ^1 



For a voltage u within the bistability range (note that 

u < and e < ) 

£/w 1 \Ef Eur , „ „, 

uth = — + — ,' <u<Uh = — A3 

67 67 1 + Q RL 67 

the electron concentration in the well is given by 
po{e^u- Ew + E^) 



Uoff = 



A + anL + 1 
po{e^u- Ew) 



(A4) 



A 



for on, intermediate and off states, respectively. The 
corresponding current- voltage characteristic is also piece- 
wise linear. 

Using (Al) and assuming for simplicity D(n) = Dq = 
const, we find an explicit solution of Eqn.(13) 

n{0 =nintexp{Pi{v)^) (A5) 
for < n < Hint and — 00 < ^ < 0; 

n-iO = iT-on - {rion - flint) exp(-/32(w)C) 

for Hint <n < Hon and < ^ < -|-oo, 

where 

/32 




Continuity of the first derivative dn/d^ at ^ = gives 
the condition 

mi \ , Q ( -.^nu-En, r. I ..e-fu-Eyj+E^ 
A A -I- URL + 1 

(A7) 

This condition (A7) implicitly determines the depen- 
dence of the front velocity v on the applied voltage u. 
Linearization of this dependence near the voltage 



Ew , 1 
Uco = h 
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that corresponds to u = leads to the explicit 
dependence (15). The obtained dependence is 
equivalent to the linearization of Eqn.(14). 
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well collector 




FIG. 1. Schematic DBRT structure. Note that the applied voltage u and the electron charge e are chosen negative. 
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2. Stationary current-voltage characteristic of the DBRT in the bistable regime. The current density J is normalized 
- eViipoE^ /h. This normalized current density is equivalent to the normalized electron concentration in the well: 
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J/ Jo = n/{poE^). The corresponding numerical values are Jo 
parameters = lOmeV, bi = 10 nm, 62 = 50nm, Tl = 0.5 meV, Fa = 0.1 meV and Fv 
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FIG. 3. Bistability range as a function of Fh/Fl- Parameters: — 10meV,6i — 10 nm, 62 ~ 50 nm and 

Fw ~ ImeV, 2meV, 3meV for curves la, 2a, 3a, respectively; = 5meV, 61 = 10 nm, 62 = 50nm and Fw ~ 1 meV, 3meV 
for curves lb, 2b, respectively; E^ — 5meV, bi — 10 nm, 62 = 30nm and Fw = ImeV, 3meV for curves Ic and 2c, 
respectively. 
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FIG. 4. Front velocity (a) and front width (b) as a function of the applied voltage it. Curves 1 and 2 correspond to Tw = 4:K 
and Tw ~ WOK, respectively. The dashed line in the panel (a) shows the prediction of the analytical formula (18) with Uco 
approximated ~ {uth + Uh)/2. Parameters as in Fig.2 : Fl = 0.5 meV, Fr = 0.1 meV, fi = 10^ cmVV • s. 
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